OFFICE  OF  NAVAL  RESEARCH 


AD-A243  200 


Grant  No.  N00014-90-J-1263 

RCT  Project  4133002 - 05 

Technical  Report  #13 


D  r  I  ( 


APPLICATION  OF  THE  MEAN  SPHERICAL  APPROXIMATION  TO 
DESCRIBE  THE  GIBBS  SOLVATION  ENERGIES  OF 
MONOVALENT  MONOATOMIC  IONS  IN  NON-AQUEOUS  SOLVENTS 


by 

L.  Blum*  and  W.R.  Fawcett* 


Prepared  for  Publication  in  the 
Journal  of  Physical  Chemistry 


*Department  of  Physics,  POB  AT,  Faculty  of  Natural  Sciences 
University  of  Puerto  Rico,  Rio  Piedras,  Puerto  Rico  00931 


**Department  of  Chemistry,  University  of  California 

Davis,  CA  95616 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government 


*This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited 


*This  statement  should  also  appear  in  Item  10  of  Document  Control 
Data  -  DD  Form  1473.  Copies  of  form  available  from  cognizant 
contract  administrator. 


91-17409 


^209  09 


'a  REPORT  SECuRiTV  CLASSiFiCAT  ON 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUM8£R(S) 

Technical  Report  #13 

6a  NAME  OF  performing  ORGANIZATION 

6d  office  S'VBOl 

Physics  Department 

University  of  Puerto  Rico 

(If  aooi'caote) 

6<.  ADDRESS  {City.  Stare,  and  ZIP  Code) 

Rio  Piedras 

Puerto  Rico  00931-3343 

3a.  NAME  OF  FUNDING 'SPONSORING 

Sb  OFF  CE  Symbol 

ORGANIZATION  Chemistry 

(If  applicable) 

Office  of  Naval  Research 

Code  472 

8c  ADDRESS  (Gty,  State,  and  ZiPCode) 

Arl i ngton 

Virginia  22217-5000 

REPORT  DOCUMENTATION  PAGE 


■  3  ve  VASk  \G5 


3  S'S's  ON  ,  avAilASiliTv  C? 


S  VONiTOR^NG  ORGANIZATION  REROR’  Ni..*/8£R'S) 


’a  NAME  OF  MONITORING  ORGANIZAT.ON 


”0  A00RE5S  iOry,  Sraf»  and  ZiPCodt) 


9  PROCUREMENT  .NSTRuMENT  DENTiPiCAT  ON  NUMBER 

RCT  Project  4133002— -05 


•0  SOURCE  OF  :-^NOiNG  NUMBERS 


PROGRAM 
£;.£MENT  NO 


PROJECT 

NO 


WORK  JNiT 
ACCESSION  NO 


11  Title  (include  Security  Clauificttion) 

Application  of  the  Mean  Spherical  Approximation  to  Describe  the  Gibbs  Solvation  Energies  of 
Monovalent  Monoatomic  Ions  in  Non-Aqueous  Solvents 


12  personal  AuThOR(S) 


13».  Type  OF  REPORT 

Summer 


L.  Blum  and  W.R.  Fawcett 


13b  TIME  COVERED 
FROM  _ TO _ 


’4  date  of  report  (rear.  Month,  Day)  iS  page  COunt 
11-09-91  32 


CDSAT;  COOES 


GROu 


'8  S.^9,£C"  'tP'-’S  Conr ''ue  on  ■*  necessary  and  identify  Oy  aiOCK  nu-roerj 


'9  abstract  ..Continue  on  reve'se  >(  necessary  and  identify  Oy  O’or»  P‘>'. 

The  mean  spherical  approximation  for  the  ion-dipole  mixture 
is  used  to  study  the  solvation  energies  of  monovalent  monoatomic 
ions  in  polar  solvents.  From  the  structure  of  the  analytical 
solution  it  is  inferred  that  there  are  two  parameters  that  should 
be  used  to  describe  the  solvation:  the  dielectric  constant  and  a 
polarization  parameter  which  plays  the  role  of  a  mean  field 
parameter  to  account  for  non-sphericity,  chemical  interactions,  and 
other  effects  related  to  failure  of  the  hard  sphere  model  for  the 
system.  The  model  was  successfully  fitted  to  data  for  the  Gibbs 
solvation  energy  of  five  alkali  metal  cations  and  their  halide  ions 
in  17  different  solvents  including  water.  On  the  basis  of  this 
analysis,  it  was  shown  that  a  single  value  of  the  polarization 


;0  D  STR  but  on  AvAiLuLBiL  "t  0=  ABS'RACT 
O  -NCLASS  F  ED'UNL'Mr?D  □  SAVE  AS  po 


::a  NAME  OF  responsible  NOw  OuAv. 

Dr.  Robert  Nowak 


DO  FORM  1473,  34  VAR  33  apf 


C  .S£PS 


;•  abs'Rac*  SEC^R.’y  C-ass-caton 


LLo  (Include  ^rea  Code)  22c  Dr~  Cz  S'VBOl 

(202)  696-4410  Code  472 


33  AP3  -a 
.1 


UNCLASSIFIED _ _ 

aecjwnr  cuxiimcAnow  or  tmi»  d«i«  gw***»o 


parameter  describes  the  data  for  the  cations  in  each  solvent,  a 
quite  different  value  being  appropriate  for  the  anions.  The 
polarization  parameters  were  found  to  be  linearly  correlated  with 
empirical  parameters  characterizing  the  solvents  basicity  (DN)  in 
the  case  oT  cations,  and  its  acidity  (E^)  in  the  case  of  anions. 


UNCLASSIFIED _ 

SCCU^^tTY  Cl. THIS  PACCr^«« 


Application  of  the  Mean  Spherical  Approximation  to  Describe  the  Gibbs  Solvation 
Energies  of  Monovalent  Monoatomic  Ions  in  Non>Aqueous  Solvents 


L.  Blum 

Department  of  Physics 
P.O.  Box  23343 
University  of  Puerto  Rico 
RioPiedras,  PR  00931-3343 


W.R.  Fawcett 
Department  of  Chemistry 
University  of  California 
Davis,  CA  95616 


i'wr 

>  'iKithi 

^  Justiric&tlwv. 


:  ly - - - 

'  DtstrlbuLlttu/ 
i  Avs t lab IIA Coi»e 


ni'U 


Abstract 


The  mean  spherical  approximation  for  the  ion-dipole  mixture  is  used  to  study  the  solvation 
energies  of  monovalent  monoatomic  ions  in  polar  solvents.  From  the  structure  of  the  analytical 
solution  it  is  inferred  that  there  are  two  parameters  that  should  be  used  to  describe  the  solvation:  the 
dielectric  constant  and  a  polarization  parameter  which  plays  the  role  of  a  mean  field  parameter  to 
account  for  non-sphericity,  chemical  interactions,  and  other  effects  related  to  failure  of  the  hard 
sphere  model  for  the  system.  The  model  was  successfully  fitted  to  data  for  the  Gibbs  solvation 
energy  of  five  alkali  metal  cations  and  their  halide  ions  in  17  different  solvents  including  water. 

On  the  basis  of  this  analysis,  it  was  shown  that  a  single  value  of  the  polarization  parameter 
describes  the  data  for  the  cations  in  each  solvent,  a  quite  different  value  being  appropriate  for  the 
anions.  The  polarization  parameters  were  found  to  be  linearly  correlated  with  empirical  parameters 
characterizing  the  solvents  basicity  (DN)  in  the  case  of  cations,  and  its  acidity  (Ej)  in  the  case  of 
anions. 
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Introduction 

Considerable  advances  have  been  made  over  the  last  thirty  years  in  the  statistical  mechanical 
description  of  electrolyte  solutions  in  molecular  solvents.  This  work  has  been  carried  out  using 
different  approaches:  computer  simulations, ^ *2  accurate  integral  equations,  such  as  the  hypemetted 
chain  equation^'^  and  simple  analytical  theories,  based  on  the  mean  spherical  approximation 
(MSA). 3  The  last  approach  is  attractive  for  the  chemists  examining  the  thermodynamic  properties 
of  electrolyte  solutions  not  only  because  the  model  gives  rather  simple  analytical  results  but  also 
because  it  shares  with  the  Debye  Hiickel  theory,  the  remarkable  simple  description  in  terms  of  a 
single  screening  parameter  for  any  arbitrary  mixture  of  electrolytes,  and  has  the  added  bonus  of 
satisfying  the  large  charge,  large  density  limits  of  Onsager.^-^ 

In  the  present  work  we  consider  a  mixture  of  a  salt  of  equal  size  hard  sphere  ions,  and  a 
solvent  which  is  represented  by  a  hard  sphere  of  different  size  with  a  permanent  point  dipole. 

Early  work  at  the  non-primitive  level^’9  was  restricted,  for  technical  reasons,  to  the  case  of  ions 
and  solvent  of  equal  diameter.  Tne  extension  of  the  analytical  result  to  mixmres  of  arbitrarily  sized 
ions  and  solvents  is  considerably  more  complex  14  5^  the  remarkable  fact  that  the  excess 
ionic  properties  depend  on  a  single  scaling,  Debye-like  parameter  is  still  retained  by  this 
approximation.  The  equations  for  the  most  general  case  appear  to  be  rather  complex,  but  the 
semirestricted  case  of  equal  size  ions  of  diameter  0\  and  a  different  size  dipole  of  diameter  On  is 
both  interesting  and  tractable. 

In  recent  years*^,  it  has  become  apparent  that  the  underlying  structure  of  the  MSA  consists  of 
a  number  of  scaling  parameters  which  is  equal  to  the  number  of  independent  interaction  parameters 
of  the  problem.  For  example,  for  the  ion-dipole  mixture,  there  are  only  two  interaction 
parameters,  the  charge  and  the  dipole  moment,  and  correspondingly,  there  are  only  two  scaling 
parameters,  F  for  the  charges,  and  X  for  the  dipoles. 

For  the  primitive  model  of  ionic  solutions  in  the  semirestricted  casein,  the  parameter  F  is 
determined  from  the  equation 
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(1) 


bT  i  1  +  r  <J  i 

where  the  ionic  charge  is  zje  and  number  density  pj  =  N[fV,  where  N,  is  the  number  of  ions  and  V 
is  the  volume  of  the  system.  The  dielectric  constant  of  the  pure  solvent  is  Ep,  the  temperature  is 
T,  and  Boltzmann's  constant  is  ke-  For  the  case  of  only  one  ion  with  diameter  Oj  the  screening 
parameter  F  is  related  to  the  Debye  screening  parameter 

47C  e^  2 
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by 


(1 +2ra)2  =  (1 +2ica)  (3) 

For  a  system  of  spheres  of  diameter  and  a  permanent  dipole  moment  |in,  the  MSA  result  can  be 

expressed  in  terms  of  a  single  parameter  X.  Following  Wertheim,'^  we  have 
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XHX  +  2)2 
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(4) 


where 

_  4;i  pnPn^  ... 

“  3kBT 

and  Pn  is  the  solvent  number  density.  Furthermore,  the  MSA  dielectric  constant  Ep  the  solvent 


is  given  by 


XHX  +  1)4 

16 


(6) 


As  has  been  often  done  in  the  literature,  the  parameter  X  can  be  computed  directly  from  the 
dielectric  constant  Ep  using  the  above  equation,  which  is  a  cubic.  This  parametrization  defines  an 
effective  polarization  parameter  which  mitigates  the  harsh  aproximations  that  go  into  the  hard 
dipole  model  of  the  solvent,  and  is  designed  to  take  into  account  polarizabilities,  softness  and  the 
non-sphericity  of  the  core. 
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In  this  paper  we  show  that  for  the  semirestricted  case  of  the  ion-dipole  model  of  ionic 
solutions,  an  ion-dipole  interaction  parameter  bi  is  obtained  from  a  cubic  equation which 
depends  on  F,  the  dielectric  constant  £„  and  the  combination 

(7) 

OiX, 

Although  both  En  and  are  functions  of  X,  they  play  different  roles.  When  On  0  we  recover 
the  primitive  model  MSA  equations,  and  the  solvation  thermodynamics  is  that  of  the  continuum 
model  as  described  by  the  Bom  equation.  The  single  parameter  appears  in  a  modified  form 
when  sticky  interactions  are  added  to  the  ion-solvent  potential.  These  sticky  potentials  could  be 
used  to  represent  specific  short  range  bonding  between  the  ion  and  the  solvent.  In  a  way,  can 
be  considered  as  a  mean  field  parameter,  similar  to  the  interaction  parameter  in  the  van  der  Waals 
theory  of  imperfect  gases.  When  the  size  of  the  solvent  is  finite,  then  we  get  different  expressions 
for  the  solvation  thermodynamics.70-12.13  These  will  be  the  subject  of  this  paper. 

Single  ion  Gibbs  solvation  energies  have  been  tabulated  by  Abraham  and  Liszi21.22^ 
Conway23  and  Marcus2^.25  for  a  large  number  of  systems.  These  authors  have  considered 
carefully  the  extra-thermodynamic  assumptions  necessary  to  estimate  single  ion  quantities  from 
experimentally  measured  properties  of  the  salts.  Standard  Gibbs  energies  of  hydration  are 
available  for  over  sixty,  chiefly  monovalent,  ions  and  Gibbs  energies  of  transfer  from  aqueous  to 
non-aqueous  solutions  reported  for  over  twenty  different  solvents.  As  a  result,  one  has  a  rich 
source  of  experimental  data  with  which  one  can  test  predictions  of  the  MSA  in  estimating  the  Gibbs 
energy  of  solvation  AGj(i). 

The  purpose  of  this  paper  is  to  examine  the  usefulness  of  the  MSA  in  accounting  for  the 
variation  in  AGj(i)  with  solvent  and  ion.  It  will  be  shown  that  using  only  two  parameters  Cn  and 
excellent  fitting  of  the  experimental  data  for  monoatomic,  monovalent  ions  is  obtained.  The 
parameter  can  be  related  in  a  more  fundamental  way  to  short  ranged,  highly  directional 
interactions  such  as  chemical  bonding.  The  consequences  of  the  present  work  are  discussed  with 
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respect  to  the  development  of  more  detailed  models  for  the  electrolyte  solution  and  applications  in 
non-equilibrium  situations. 

On  the  basis  of  the  MSA,  expressions  have  been  developed  for  the  ion-dipole  and  dipole- 
dipole  contributions  to  the  Gibbs  free  energy  of  solvation. An  excellent  fit  between  the 
tabulated  data  and  the  estimates  of  the  model  is  obtained  using  one  adjustable  parameter,  the 
polarization  parameter  X.  It  is  shown  that  the  values  of  X  which  bring  about  agreement  between 
theory  and  experiment  are  roughly  correlated  with  the  solvent's  basicity  in  the  case  of  cations  or 
with  the  solvent's  acidity  in  the  case  of  anions.  These  results  can  be  interpreted  within  the  MSA 
as  indicating  that  specific  short  range  interactions  (chemical  association  effects)  are  impOTtant  The 
present  paper  describes  a  method  of  including  these  "sticky"  interactions  within  the  context  of  the 
MSA  theory  on  the  basis  of  an  analysis  of  data  for  the  Gibbs  solvation  energies  of  monovalent 
monoatomic  ions  in  17  solvents,  both  protic  and  aprotic. 
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Theory 


We  summarize  the  results  of  previous  work.'°  W’e  use  the  invariant  expansion  fonnalisiu. 
in  which  the  total  pair  correlation  h{l2)  is  expanded  in  terms  of  rotational  invariants 

h{12)  =  h^^iryi)  +  +  ^"^(^12)^'"  (8) 

where  h"’"^{ri2)  is  the  coefficient  of  the  invariant  expansion,  which  depend  only  on  the 
distance  ri2  between  spheres  1  and  2.  The  rotational  invariants  depend  only  on  the 

mutual  orientations  of  the  molecules.  For  the  present  case  the  relevant  correlation  functions 


•  lon-ion: 


/,..(r)  =  (l/2)te{r)-^°oe{r) 


•  ion-dipole: 


/»,„(;■)  =  (1/2)  [//!'„'('•) 


•  dipole-dipole; 

h„n{r)  =  -v''3/?,\;f(r)//,  •  i)> 

+  -/'i  -  Pi]  (H) 

where  fi  is  the  unit  vector  in  the  direction  of  //.  The  solution  of  the  MS.A  is  given  in 
terms  of  the  'energy'  parameters 


•  lon-ion: 


•  ion-dipole: 


bo  =  2xp,  / 
Jo 


ilrh„{r)r 


•  dipole-dipole: 


drh,„{r) 
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which,  as  will  be  shown  below  are  pruj  rional  to  the  ioii-ion.  ion-dipole  and  dipole-dipole 
excess  internal  energy."  In  the  MSA  they  are  functions  of  the  ion  charge  and  the  solvent 
dipole  moment,  through  the  parameters 


hsT 


(15) 


and 

a  _  d7rp„p^ 

3/.sr 

These  parameters  are  required  to  satisfy  the  following  etjuation-s'* 


(16) 


c/f  -f  o  \  —  (Iq 


~  b  ^  fii  + 

Z< 


where  w'e  have  used  the  definitions 
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+  "r"i 


U2Dy 


Df^ 


h\rr, 

■^(1  4-  ho)  -f 


Urr, 


h]  , 
A  =  —  A  d* 


and 


Also 


A^di- 


hldoO’n 

2A<t, 


1  A  ho  dpcr^ 
2  ^  6(7, 


4  =  1-  h,/6 

•d}  —  \ 


.Vi  = 


(1  Ah,/12)^ 

As  has  been  previ  ,ly  shown'*’  Eq.(18)  can  be  replaced  by 

Oihirr„ 


2/7, 


A  n  1  d»  —  doAi/i 


(17) 

(18) 

(19) 


(- 

(21) 

(22) 

(23) 

(24) 

(25) 
(2o, 

(27) 

(28) 
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which  together  with  Eqs.(li)  and  (19)  form  a  system  of  thre*>  equations  for  the  three  un¬ 
knowns  bo,  bx,b2.  However  these  equations  are  stiU  complex.  .A  much  simpler  s.  r  of  e.piations 
IS  obtained  when  we  use  the  proper  scaling  lengths'”  We  define  T  through  the  relation 


U-1- 


^  +  Ter,  4  1*  1  -H  Tct,  3rr, 


Ter,  = 


and  A  from 


- U.:-  - 


After  some  long  but  straightforward  algebra  we  get 


dx  =  ^rcr,(l  +  Trr, 


(1 2  =  (li 


2^C. 


where  we  have  defined  the  convenient  parameter  c, 


1  +  Ter,  -t- 

tTr, 

1  -  r^,  +  1 


v=  1 


Furthermore  we  get  for  Eq.(24) 

A  =  (I2  -  </,] 

and  the  MSA  equation  28  reduces  to 


'>r(  A  -  2)* 


\  f  - 

12.1J  \  1  +  r^, 


2rcT,(i  +  -(p-  = 


V 

together  witu  the  subsidiary  cubic  equation  for  the  ion-dipole  parameter  fi, 

P'[c/  +  D  -  Ij  =  [C.  -  (p  -  1)C  j 
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(39) 


which  can  be  written  in  a  more  convenient  form  using 


„  _  _  hi(\^2)^ 

-  TTT  - 


36 


Then.  Eq.(37)  reads 


2r-7,(i  +  T/t.)  [1  -  ei(^./c,)]  = 


</n(l  *  B,y 


V 


and  Eq.(38)  is  a  cubic  equation  in  B\ 


(40) 


( 1  +  Si )‘  |c,‘  +  Si  j  —  6n  is.  “  Bx^,,Y  (41 ) 

which  is  our  most  interesting  new  result  for  the  ion-dipole  MSA.  This  equation  is  considerably 
simpler  than  the  ones  derived  pre^^ously.  since  the  ion-dipole  coupling  parameter  6],  or  Si. 
which  is  the  responsible  for  the  deviations  of  the  molecular  solvent  model  from  the  primitive, 
continuum  dielectric  model  appears  to  be  a  function  of  the  dielectric  constant  and  the 
ionic  screening  parameter  F.  If  we  let  the  diameter  of  the  solvent  shrink  to  zero,  then  from 
Eq-(~)  SI  — ’  0  .  and  Eq.(41)  becomes 

(l+Si)'S,  =0  (42) 

which  has  only  one  real  solution  ^i  =  0.  There  is  no  (onpling  between  the  ions  and  the 
dipoles,  and  Eq.(40)  reduces  to 

2ra-,(l  -t-  F'T,)  =  (43) 

V 

which  should  be  compared  to  the  primitive  model  MS.\  equation  Eq.3  in  the  form 

2ro-(l  +  Ttr)  =  Krr  I  44 ) 

From  the  results  of  the  sticky  model'^  we  know  that  when  short  ranged  interactions  are 
present,  they  change  .f,  and  the  dielectric  constant  e  in  a  different  way.  This  means  that  if 
we  are  thinking  of  our  model  as  a  simple  representation  of  a  real  solvent,  the  two  parameters 
should  be  treated  independently.  The  connecting  equation  (6)  between  the  parameters  \ 
and  the  dielectric  constant  en  is  strictly  a  consequence  of  the  oversimplified  spherical  dipole 
representation  of  the  solvent.  more  general  mean  field  (  in  the  van  der  Waals  equation 
sense)  model  should  treat  both  of  these  parameters  independently,  and  consequently,  assign 
their  frequency  behaviour  accordingly.  This  is  what  is  being  done  in  the  next  sei  tion  to 
study  the  solvation  thermodynamics  for  a  number  of  solvents. 
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A  second  limit  of  interest  is  the  infinite  dilution  limit.  This  was  first  done  by  Blura'*^ 
who  derived  the  limiting  expression  for  the  excess  internal  energy  parameter  hi.  It  was  later 
discussed  in  much  more  detail  by  Chan  and  colaborators,'’’  and  in  a  more  comprehensive 
way  by  Garisto  et  al.*°  From  Eq.(41) 


Results  and  Discussion 

In  order  to  test  the  application  of  the  MSA  expression  for  the  Gibbs  energy  of  ion  solvation, 
one  must  choose  experimental  data  which  reflect  the  properties  of  single  ions  under  standaid 
conditions.  Extensive  lists  of  thermodynamic  properties  in  water,  and  for  transfer  from  water  to 
various  non-aqueous  solvents  have  been  compiled  by  Abraham21.22^  Conway^^,  and  Marcus.2'*.25 
Separation  of  the  thermodynamic  properties  of  electrolytes  into  components  for  the  cation  and 
anion  is  a  well  known  problem  in  solution  electrochemistry.  It  has  been  discussed  in  detail  by 
Conway23.27^  ^ho  has  assessed  the  relative  merits  of  the  commonly  used  extrathermodynamic 
assumptions  especially  for  aqueous  solutions.  We  avoid  consideration  here  of  the  best  source  of 
Gibbs  energies  of  solvation  for  single  ions,  and  choose  the  values  given  by  Conway^^  simply 
because  these  values  give  the  best  fit  with  the  present  theory. 

The  standard  Gibbs  energy  of  solvation  of  ion,  i,  in  solvent  n,  AG°  (i),  was  estimated  from 
the  same  quantity  measured  in  water,  AG^  (i),  using  tabulated  values  of  the  Gibbs  energy  of 
transfer,  AG^  (i)  from  water  to  the  second  solvenv.  Thus,  the  relationship  between  these  quantities 
is 

AGj(i)  =  AG^(i)  +  AG°  (i)  (54) 

An  extrathermodynamic  assumption  is  also  involved  in  obtaining  single  ion  Gibbs  energies  of 
transfer.  The  commonly  accepted  assumption^*'^^  is  that  the  Gibbs  energy  of  transfer  for  the  salt 
tetraphenylarsonium  tetraphenylborate  may  be  divided  into  equal  contributions  for  cation  and  anion 
(TATB  assumption). 

Separate  fits  of  theory  to  thermodynamic  data  were  performed  for  cations  and  anions.  The 
cations  considered  were  the  five  alkali  metal  ions  for  which  extensive  data  are  available,  namely, 
Li"*",  Na+,  K+,  Rb'*’  and  Cs'*'.  The  monoatomic  anions  for  which  equally  numerous  data  can  be 
found  are  Ck,  Br,  and  k.  The  fluoride  ion  was  not  considered  because  its  salts  are  not  very 
soluble  in  non-aqueous  solvents  and  very  few  thermodynamic  data  have  been  tabulated  for  this 
ion. 25  The  ionic  radii  needed  to  estimate  the  Gibbs  energies  were  assumed  to  be  the  Pauling 
values.  The  values  of  AG^  (i)  used  together  with  these  radii  are  given  in  Table  1 . 
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As  descnbed  above,  the  thennodynamic  properties  of  a  single  ion,  i,  may  be  estimated  in  a 
given  solvent  using  the  MSA  with  knowledge  of  two  parameters  and  4i-  In  the  following 
analyses,  En  has  been  always  set  equal  to  the  solvent’s  static  dielectric  constant.  The  parameter 
may  be  resolved  into  two  components,  Oj,  the  diameter  of  the  ion,  and  X/On,  which  we  will  show 
depends  only  on  the  solvent  and  the  sign  of  the  charge  on  the  ion.  An  initial  estimate  of  X/On  may 
be  obtained  by  calculating  X  on  the  basis  of  eq.  (6)  and  using  values  of  On  estimated  on  the  basis  of 
gas  solubilities.28  This  estimate  gives  values  of  AGn  (i)  which  are  close  to  those  estimated  from 
experimental  data  for  the  monoatomic  monovalent  ions  considered  here.  In  order  to  obtain  a  better 
fit  between  experimental  data  and  theory,  A/On  was  varied  in  each  solvent,  cations  being 
considered  separately  from  anions.  The  estimate  of  AGn  (i)  was  based  on  the  ion-dipole 
contribution  only  (eq.  52),  the  dipole-dipole  contribution  being  ignored;  that  is,  we  assume  that 
Gdd  is  seriously  overestimated by  eq.  53.  Thus,  our  estimates  of  AGn,th  (i)  are  based  on  the 
equation 

AGn.th  (i)  =  -  (55) 

where  Nq  is  Avogadro's  number  and  is  given  by  eq.(7).  The  criterion  used  to  determine  the 
best  fit  was  a  minimum  in  the  standard  deviation  defined  as 

s.d.  =  (X  [AG°  g^(i)-  AG°  jh(i)P/(n-l))’/2  (56)  • 

i  ‘ 

Values  of  X/Cn  were  determined  in  17  different  solvents  including  water  on  the  basis  of 
reported  data  for  the  five  alkali  metal  cations.  The  fit  of  theory  to  experiment  was  excellent  with 
an  average  standard  deviation  of  7.6  kJ  mol'*.  The  solvents  considered  were  restricted  to  those 
which  are  reasonably  polar,  that  is,  with  dielectric  constants  greater  than  20.  The  resulting  values 
of  X/On  are  summarized  in  Table  2.  This  parameter  is  close  to  six  for  most  solvents  and  usually 
greater  than  the  estimate  based  on  values  of  X  calculated  using  eq.  (6),  which  are  between  2  and 
2.5  and  values  of  On  obtained  from  gas  solubilities^*,  which  are  close  to  0.5  nm.  Thus,  values  of 
X/(5n  on  the  basis  of  the  simple  model  are  typically  between  4  and  5.  The  present  result  that  the 
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values  obtained  on  the  basis  of  a  fit  to  experimental  data  are  higher  is  attributed  to  the  fact  that  the 
simple  model  does  not  account  for  the  specific  chemical  interactions  between  the  cation  and  the 
solvent  dipole  which  contribute  to  its  solvation. 

In  order  to  illustrate  the  quality  of  the  fit  of  the  present  model  to  experiment,  values  of  the 
estimated  free  energy  of  transfer  are  plotted  against  the  experimental  quantity  for  three  of  the 
cations  considered  in  Figs.  1-3.  It  should  be  noted  that  this  method  of  presenting  our  results  is  the 
most  demanding  since  AG^  constitutes  only  a  fraction  of  the  total  free  energy  of  solvation.  In  this 
regard,  the  estimated  value  of  the  Gibbs  energy  of  transfer  is  defined  by  the  equation 

=  ^^Ith  -  ^G;,ex  (57) 

where  AG°  is  the  estimated  value  of  the  Gibbs  energy  of  solvation  in  solvent  n.  In  the  case  of 
Na'*'  (Fig.  1)  the  fit  is  excellent,  the  difference  between  the  estimated  and  experimental  values  of 
AG^  being  only  a  few  kJ  mol'^  Similar  results  are  obtained  for  the  K"*"  ion,  but  in  this  case  the 
estimated  values  tend  to  be  lower  than  the  experimental  ones  by  a  few  kJ  mol'^  Results  for  the 
Li'*’  system  shown  in  Fig.  2  show  greater  deviations  from  experiment,  the  deviations  in  this  case 
being  usually  positive.  The  poorest  fit  was  obtained  for  Cs'*'  (Fig.  3)  where  negative  deviations 
approaching  20  kJ  mol'i  are  seen  for  some  solvents.  Some  of  the  discrepancy  between  theory  and 
experiment  may  be  due  to  errors  in  the  extrathermodynamic  assumptions  used  to  extract  single  ion 
solvation  energies.  However,  the  main  reason  for  the  variation  in  the  quality  of  fit  with  the  nature 
of  the  cation  is  the  fact  that  the  present  analysis  ignores  the  contribution  to  the  Gibbs  energy  from 
dipole-dipole  interactions.  These  will  increase  in  importance  with  decrease  in  ion  size  because  the 
strength  of  the  ion-solvent  bonds  for  small  ions  like  Li'*’.  The  other  parameter  which  could  be 
varied  to  improve  the  fit  is  the  ionic  radius,  which  has  been  arbitrarily  set  equal  to  the  Pauling 
value.  This  assumption  is  considered  funher  below.  When  one  keeps  in  mind  that  the  relative 
error  for  the  total  solvation  energy  is  much  smaller  than  that  for  the  Gibbs  transfer  energy,  it  is 
clear  that  the  MSA  model  gives  a  very  good  description  of  the  experimental  solvation  energies  for 
the  alkali  metal  cations  using  one  adjustable  parameter  for  each  solvent,  namely,  VcTn. 
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A  similar  procedure  was  used  to  fit  the  Gibbs  solvation  energy  data  for  the  halide  ions.  Cl', 
Br,  and  I'.  In  this  case,  the  average  standard  deviation  of  the  fit  between  theory  and  experiment 
was  8.3  kJ  mol-'.  For  anionic  solvation,  the  estimated  values  of  the  parameter  X/On  are  very 
different  from  those  obtained  for  the  cations  (see  Table  2).  For  the  aprotic  solvents,  this  quantity  is 
usually  in  the  range  from  10  to  12  nm'i  whereas  for  protic  solvents  it  is  between  17  and  24  nm'*. 

It  follows  that  the  MSA  estimates  of  solvation  free  energies  for  anions  are  much  closer  to  the  Bom 
estimate  than  those  for  cations.  As  Von  becomes  larger,  the  quantity  is  smaller,  and  the  MSA 
factor,  (1  +  ^i)'i  in  eq.  (55)  becomes  less  imponant  This  result  is  not  unexpected  because  it  is 
well  known  that  typical  polar  solvents  can  effectively  solvate  cations  through  the  lone  electron  pairs 
on  the  electro-negative  group  in  the  solvent  molecule,  whereas  anion  solvation  is  much  weaker 
unless  the  solvent  is  protic  and  can  stabilize  the  anion  through  hydrogen  bonding. 

A  plot  of  AG^  against  AG^  is  shown  for  the  Cl'  ion  data  in  Figure  4.  Agreement 
between  the  theoretical  estimate  and  experiment  is  quite  good  except  for  higher  values  of  AG° 
where  positive  deviations  are  found.  In  the  case  of  Br  ion,  the  variation  in  AG^  follows  that 
in  perfectly  except  that  the  theoretical  estimates  are  high  by  about  8  kJ  moF’  (Fig. 

5).  These  results  suggest  that  the  value  of  AG^  (Br)  may  be  in  error  by  8  kJ  mol'F  Finally,  the 
poorest  fit  between  theory  and  experiment  is  found  for  the  I'  data  where  negative  deviations  in 
AG^^jjj  are  seen  for  lower  values  of  AG°^j^.  However,  when  one  considers  that  the  fitting 
was  achieved  with  one  adjustable  parameter  and  that  the  relative  error  in  AG°  (i)  is  much  smaller, 
the  results  obtained  are  remarkably  good. 

As  pointed  out  above,  the  values  of  the  parameter  X/On  obtained  in  fitting  the  data  for  the 
Gibbs  energy  of  transfer  differ  significantly  from  those  predicted  on  the  basis  of  the  simple  model 
(eq.  6),  and  differ  between  cations  and  anions.  Since  one  expects  that  this  parameter  contains 
some  measure  of  the  chemical  ability  of  the  solvent  to  solvate  ions,  it  is  interesting  to  see  whether 
there  is  any  connection  between  the  parameters  obtained  here  and  empirical  parameters  used  to 
measure  the  solvent's  ability  to  solvate  cations  and  anions.  In  the  case  of  cations,  a  well  known 
parameter  measuring  the  solvents  basicity  or  cation  solvating  ability  is  the  donor  number  (DN) 
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introduwed  by  Gutinann.29  On  the  other  hand,  the  solvent’s  anion  solvating  ability  or  acidity  is 
empirically  measured  by  the  Dimroth-Reichardt  parameter  Et-^  It  was  shown  earlier^'-^^  that 
Gibbs'  solvation  energies  and  other  themx)dynamic  parameters  for  ions  in  electrolyte  solutions  are 
linearly  related  to  these  parameters  by  the  equation 

Q  =  Qo  +  aEj  +  pDN  (58) 

where  Q  is  the  thermodynamic  quantity  in  question,  Qq,  the  solvent  independent  value,  a,  a 
quantity  measuring  response  to  solvent  acidity,  and  j3,  that  to  solvent  basicity.  Since  effects  for 
cations  and  anions  have  been  separated  in  the  present  analysis,  it  is  reasonable  to  expect  that  only 
one  of  these  parameters  is  important  in  explaining  the  variation  in  XyOn  with  solvent  nature. 

A  plot  of  Von  for  cations  against  the  solvent's  donor  number  is  shown  in  Fig.  7.  A  very 
good  correlation  is  obtained  using  data  for  the  1 1  aprotic  solvents  involved  in  the  present  study 
with  a  correlation  coefficient  of  0.928.  It  should  be  noted  that  data  are  available  for  solvents  with 
a  wide  range  of  donor  properties  from  very  poor  Lewis  bases  like  nitromethane  and  nitrobenzene 
to  strong  Lewis  bases  such  as  hexamethylphosphoramide.  The  present  results  provide  str 
confirmation  that  the  values  of  X/On  obtained  here  by  fitting  data  for  the  alkali  metal  cations  no 
account  for  the  specific  chemical  interaction  between  the  electro-negative  group  on  the  solvent 
molecule  and  the  cation.  As  is  well  known2^>29,  the  donor  number  is  not  well  defined  for  protic 
solvents  so  that  these  data  could  not  be  included  in  the  present  correlation.  However,  in  the  case 
of  water,  the  donor  number  is  often  assumed  to  be  18.  The  present  result  for  water  fits  well  with 
those  for  other  solvents  using  this  value. 

A  plot  of  the  values  of  X/<Jn  obtained  by  fitting  the  data  for  the  halide  ions  against  the 
solvent’s  acidity  parameter,  Ey  is  shown  in  Fig.  8.  In  this  case  results  for  16  solvents  are 
available,  and  a  very  good  correlation  is  obtained  (r  =  0.944).  It  is  apparent  that  the  properties  of 
the  correlation  are  determined  to  a  great  extent  by  the  protic  solvents  which  have  both  high  values 
of  Ej  and  VOn.  The  values  of  Ej  for  the  aprotic  solvents  are  close  together,  the  corresponding 
values  of  X/On  falling  in  the  range  from  9  to  13  nm'L  A  much  larger  variation  in  VOn  is  found 
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for  anionic  solvation  and  it  clearly  is  related  to  the  solvent's  acidity  as  estimated  by  the  parameter 
Et. 

Having  shown  that  reasonable  fits  between  theory  and  experiment  are  possible  using  one 
adjustable  parameter,  and  that  this  parameter  is  related  to  well  known  empirical  parameters 
describing  solvent  basicity  and  acidity,  it  is  interesting  to  test  the  ability  of  the  present  model  to 
predict  Gibbs  solvation  energies  for  other  ions.  We  have  carried  out  this  test  in  two  ways, 
namely,  by  examining  data  for  monoatomic  divalent  ions,  and  for  non-spherical  monovalent  ions. 

A  plot  of  the  Gibbs  energy  of  solvation  for  the  alkaline  earth  metal  ions  together  with  data  for 
Zn2+,  Cd2+  and  Pb2+  ions  is  shown  in  Fig.  9.  These  estimates  were  also  made  assuming  the 
Pauling  radii  for  the  ions.  It  is  clear  that  the  present  theory  gives  an  excellent  estimate  of  AG^  (i) 
in  aqueous  solutions  for  Mg2+,  Ca2+,  Sr2+  and  Ba2+  but  that  it  underestimates  this  quantity  for  the 
other  ions  by  approximately  200  kJ  mol'f  The  failure  for  ions  such  as  Cd2+  is  attributed  to  the 
fact  that  the  bonds  between  the  water  molecules  and  ions  involve  d  electrons.  Thus,  the  value  of 
VOn  derived  here  is  not  appropriate  for  transition  metal  cations.  Gibbs  energies  of  transfer  from 
water  to  non-aqueous  systems  are  available  for  Ba2+  in  some  solvents.25  Values  of  AG^ 
estimated  for  this  ion  using  the  present  model  agree  quite  well  with  those  obtained  from  the 
experimental  data.24.25 

An  attempt  was  also  made  to  apply  the  present  model  to  data  for  non-spherical  ions.  In  the 
case  of  C104‘,  a  tetrahedral  ion,  the  estimate  of  AG^  in  water  is  -262  kJ  mob'  assuming  an  ionic 
radius  of  0.241  nm  .33  The  estimated  value  from  experimental  data  is  -214  kJ  mob'.  24  Thus,  the 
present  model  significantly  overestimates  the  Gibbs  solvation  energy  for  small  tetrahedral  ions. 
This  suggests  that  the  mechanism  of  ionic  solvation  is  significantly  different  from  that  for 
monoatomic  anions,  the  resulting  value  of  X/On  being  smaller.  In  the  case  of  tetrahedral  cations, 
the  situation  is  somewhat  better.  The  present  model  gives  an  estimate  of  AG°  equal  to  -297  kJ 
mob'  for  NH4+  assuming  a  radius  of  0.148  nm,  whereas  the  corresponding  experimental  estimate 
is  -292  kJ  mob'. 24  in  the  case  of  (CH3)4N'*’  and  (C2H5)4N^  the  corresponding  results  are 
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-189  and  -163  kJ  mol‘*  from  theory,  and  -200  and  -172  kJ  mol"^  from  experiment. 

Comparison  of  theoretical  and  experimental  results  suffers  for  these  systems  since  the  experimental 
quantities  have  been  estimated  by  several  authors  using  different  extrathermodynamic  assumptions. 
However,  it  would  seem  that  the  estimates  from  the  present  model  are  reasonably  close  to 
experimental  values  for  aqueous  solutions. 

It  is  interesting  to  speculate  on  the  significance  of  the  TATB  assumption  with  respect  to  the 
present  results.  On  the  basis  of  the  results  for  monoatomic  cations  and  anions,  the  values  of  X/rtn 
are  so  different  that  one  would  not  be  able  to  choose  large  spherical  ions  for  which  AG°  would  be 
equal  for  a  cation  and  anion  of  the  same  size.  However,  the  brief  examination  of  data  for 
tetrahedral  ions  made  here  suggests  that  X/Gn  appropriate  for  tetrahedral  cations  and  anions  may  be 
much  closer  to  one  another  than  for  monoatomic  ions.  In  this  case,  the  MSA  model  would  predict 
that  values  of  AG°  for  these  ions  are  equal,  and  the  TATB  assumption  valid. 

Some  comment  should  be  made  here  about  the  values  chosen  for  the  radii  of  the  ions.  It  is 
clear  from  previous  discussion  in  the  literature23.33^  and  also  from  our  work,  that  the  choice  of 
ionic  radius  is  very  important  in  determining  the  magnitude  of  the  estimate  of  AG^  (i),  and,  in  the 
present  case,  the  value  of  the  MSA  parameter  A/Op.  In  the  present  study,  Pauling  radii  were  used 
for  monoatomic  ions.  Otherwise,  the  radii  tabulated  by  Marcus^'^  were  chosen.  In  more  recent 
work,  Marcus33  has  recommended  using  revised  radii  obtained  from  crystallographic  data,  mainly 
for  ions  with  a  coordination  number  of  six,  as  estimated  by  Shannon  and  Prewitt.^*^  The  same 
author  has  also  tabulated33  ionic  radii  estimated  from  X-ray  data  for  aqueous  solutions.  However, 
the  best  fit  to  the  present  data  was  obtained  for  the  alkali  metal  cations  using  the  original  Pauling 
radii.  Variation  in  the  thermodynamic  projjenies  of  single  ions  with  ionic  size  has  been  used  as  a 
method  of  assessing  the  extrathermodynamic  assumptions  made  to  separate  cationic  and  anionic 
contributions  so  that  the  arguments  used  to  select  either  the  ionic  radii  or  the  extrathermodynamic 
assumptions  are  not  independent  of  one  another.  However,  on  the  basis  of  the  single  ion  data 
currently  used  in  the  literature^'^’^^  the  original  Pauling  radii  are  preferrable. 
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In  summary,  the  approach  used  here  in  which  the  parameter  X/On  is  estimated  from  the  large 
body  of  experimental  data  available  for  monoatomic  ions  has  proven  to  be  remarkably  successful. 
It  shows  clearly  that  the  chemical  interaction  between  the  ion  and  solvent  dipole  needs  to  be 
considered  in  developing  a  reasonable  statistical  mechanical  model  for  electrolyte  solutions.  In 
future  papers,  we  will  consider  other  thermodynamic  properties  of  these  systems. 
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Table  1.  Pauling  radii  and  standard  Gibbs  energy  of  solvation  of  monovalent  ions  in  water  at 
298°  K23 


Ion 

Pauling  Radius 

ri,  nm 

Li+ 

0.060 

Na+ 

0.095 

K+ 

0.133 

Rb+ 

0.148 

Cs+ 

0.169 

Cl- 

0.181 

Br 

0.195 

I- 

0.216 

Standard  Gibbs  energy  of  hydration 
AG^  (i),  y  mol-i 
-487.5 
-387.6 
-314.0 
-292.7 
-260.5 

-339.7 

-325.9 

-279.5 
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Table  2.  '^dlues  of  the  MSA  parameter  XyOn  for  monovalent  monoatomic  cations  and  anions 
obtained  by  fitting  experimental  data  to  eo.  55  in  17  different  solvents 
Solvent  Vcn,  nm-i 


Cations 

Anions 

1. 

Water  (w) 

6.06 

24.00 

2. 

Methanol  (MeOH) 

6.04 

20.48 

3. 

Ethanol  (EtOH) 

5.96 

18.58 

4. 

1 -Propanol  (PrOH) 

6.02 

17.36 

5. 

Formamide  (F) 

6.26 

17.06 

6. 

N-methylformamide  (NMF) 

6.35 

- 

7. 

Acetone  (AC) 

6.27 

11.01 

8. 

Acetonitrile  (AN) 

5.74 

12.20 

9. 

Dimethylacetamide  (DMA) 

6.75 

10.04 

10. 

Dimethylformamide  (DMF) 

6.60 

11.10 

11. 

Dimethylsulphoxide  (DMSO) 

6.64 

12.77 

12. 

Hexamethylphosphoramide  (HMPA) 

6.89 

9.61 

13. 

Nitrobenzene  (NB) 

5.31 

13.18 

14. 

Nitromethane  (NM) 

5.33 

12.76 

15. 

N-methylpyrrolidinone  (NMP) 

7.01 

11.04 

16. 

Propylene  carbonate  (PC) 

5.71 

11.93 

17. 

Tetramethylene  sulphone  (TMS) 

6.24 

11.01 
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Figure 
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Figure 

Figure 

Figure 

Figure 

Figure 


Legends 

1.  Plot  of  the  theoretical  estimate  of  the  Gibbs  energy  of  transfer  of  Na"*"  ion  from 
water  to  a  non-aqueous  solvent  against  the  corresponding  value  from  experiment. 
The  integers  indicate  the  solvent  according  to  the  list  given  in  Table  2.  The  straight 
line  is  drawn  with  unit  slope. 

2.  As  in  Fig.  1  but  for  the  Li"^  ion. 

3.  As  in  Fig.  1  but  for  the  Cs'*’  ion. 

4.  As  in  Fig.  1  but  for  the  Cl*  ion. 

5.  As  in  Fig.  1  but  for  the  Br*  ion. 

6.  As  in  Fig.  1  but  for  the  r  ion. 

7.  Plot  of  the  MSA  parameter  X/Og  estimated  for  monovalent  cations  in  various  non- 
aqueous  solvents  as  indicated  against  the  solvent's  donor  number  DN. 

8.  Plot  of  the  MSA  parameter  X/Og  estimated  for  monovalent  anions  in  various  non- 
aqueous  solvents  as  indicated  against  the  solvent’s  acidity  parameter,  Ex- 

9.  Plot  of  the  Gibbs'  solvation  energy  for  various  divalent  cations  in  water  against  the 
ionic  radius.  The  curve  shows  the  estimate  of  the  same  quantity  according  to  the 
MSA  model  using  X/Og equal  to  6.06  nm*’. 
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